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CONSPECTUS 

Interest in electrochemistry is surging, driven by new applications in energy conversion, water 
treatment, materials processing, and biotechnology. As geometries shrink to the nanoscale, the 
rate-limiting step is often ion intercalation (i.e. reversible insertion) into a host solid for transport 
or storage. For example, oxygen intercalates into a ceramic electrolytes in solid oxide fuel cells, 
and lithium intercalates into carbon or metal oxide nanoparticles in Li-ion batteries. The standard 
phenomenological model for electrode kinetics is the Butler- Volmer equation, which fits the current- 
voltage relation in many situations and can be justified (in certain limits) by the Marcus theory of 
charge transfer. Existing theories, however, provide little guidance as to the form of the exchange- 
current prefactor to account for configurational entropy, elastic stress, phase transformations, and 
other non-idealities arising in ion intercalation. 

These challenges are exemplified by the high-rate cathode material, Li^FePCU (LFP), which has 
a strong tendency to phase separate into Li-rich and Li-poor phases that is believed to limit its 
performance. Phase separation was originally thought to occur as an isotropic "shrinking core" 
in each particle, but experiments later revealed striped phase boundaries along the active facet. 
Meanwhile, dramatic rate enhancement was attained with LFP nanoparticles, and classical battery 
models could not predict the role of phase separation. This Account describes the development of 
a theory that answers this question via a variational formulation of Faradaic reaction kinetics for 
ionic solids and concentrated solutions. 

The theory is based on non-equilibrium thermodynamics, consistent with Cahn-Hilliard phase- 
field models for the solid host. Butler- Volmer and Marcus kinetics are reformulated for concentrated 
solutions using activity coefficients. The theory is applied to lithium insertion in transition metal 
oxides. For phase-separating solids, such as LFP, the intercalation rate is enhanced by concentra- 
tion gradients and elastic coherency strain. This causes exposed phase boundaries to propagate as 
"intercalation waves" at low current. Above a small critical current, homogeneous reactions are 
favored, which helps to explain the high rate capability of LFP nanoparticles. The theory also 
predicts similar phenomena in porous electrodes with phase-separating nanoparticles. Narrow re- 
action fronts with mosaic instabilities at low currents become broadened and limited by electrolyte 
diffusion at high currents. 



I. INTRODUCTION 

The standard model of electrode kinetics is the Butler- 
Volmer (BV) equationp] E], 

I = Iq ^ e - a cnerj/k B T _ e a a ner)/k B T^j ^) 

which relates the total current / due to the Faradaic reac- 
tion, O + ne~ — ^ R, which converts the oxidized state O 
to the reduced state R and consumes n electrons, to the 
exchange current Iq and overpotential 77 = A</> — A</> eg , 
where Acj) is the interfacial voltage whose Nernst equi- 
librium value is A(j) eq . For a simple single-step charge- 
transfer reaction, the anodic and cathodic charge-transfer 
coefficients a a and a c should satisfy a a = 1 — a and 
a c = a with a symmetry factor, < a < 1, although 
the same form is often fitted to more complicated exper- 
imental voltage dependence. The microscopic theory of 
charge transfer [3], pioneered by R. A. Marcus and hon- 
ored by the Nobel Prize in Chemistry [4j, provides justi- 
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fication for the BV equation and a means to estimate its 
parameters in dilute solutions pQ. 

Since ion intercalation occurs concentrated solid solu- 
tions, the reaction rate is affected by packing entropy, 
elastic strain, and other non-idealities, which are ne- 
glected in classical theories. In Li-ion battery simula- 
tions [2j [5HZ] ? only excluded volume has been considered, 
by empirically modifying the exchange current 

Io(c) = k{c max - c) a ^c a ' (2) 

where c max is the maximum lithium concentration in the 
solid. Moreover, the Nernst potential Acj) eq (c) is fitted 
to the open circuit voltage, and the chemical diffusivity 
D chem (c) in the solid is adjusted to fit discharge curves. 
In a thermodynamically consistent theory, however, these 
properties must all derive from the local lithium activity 
at the nanoscale. 

Advances in Li-ion batteries increasingly rely on the 
unusual kinetics of nanoparticles, as best illustrated by 
lithium iron phosphate, Li x FeP04 (LFP). In the seminal 
LFP paper, Padhi et al. [8] concluded that "the material 
is very good for low-power applications" but "at higher 
current densities there is a reversible decrease in capacity 
that... is associated with the movement of a two-phase 
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interface" (between Li-rich and Li-poor phases). Over 
the next decade, LFP ironically became a popular high- 
power material [9HTT], but only in nanoparticle form. Un- 
derstanding the role of phase separation in this reversal 
of fortune turned out to be a major scientific challenge. 

It is now understood that phase separation is strongly 
suppressed in LFP nanoparticles, to some extent in equi- 
librium [T2HT4] , but especially under applied current [14]- 
[17], since reaction limitation [18], anisotropic lithium 
transport p~9j [20] , elastic coherency strain [14] [2TH23] , 
and surface energy [15] [24] [25] are enhanced. At low 
currents, however, anisotropic nucleation and growth 
can also occur in nanoparticles [14] [15] [18] [26] , leading 
to stochastic multi-particle intercalation [27-30 . These 
phenomena could not be described by traditional battery 
models [2] [5] , which assume a spherical "shrinking core" 
phase boundary [6] [7] . 

This Account summarizes my struggle to develop 
a theory of intercalation kinetics in phase-separating 
electrodes [T5HT5] fT8l [23] [30H55] by combining charge- 
transfer theory [3 with concepts from statistical 
physics [36] [37] and non-equilibrium thermodynam- 
ics [38j |39] . The work began when my postdoc, G. Singh, 
found the paper of Chen et al. [40] revealing striped phase 
boundaries in LFP, looking nothing like a shrinking core 
and suggesting phase boundary motion perpendicular to 
the lithium flux, which could not be described by stan- 
dard electrochemical kinetics. G. Ceder suggested adapt- 
ing the Cahn-Hilliard (CH) model for bulk phase sepa- 
ration [41], but it took several years to achieve a con- 
sistent theory. Our initial model predicted intercalation 
waves [18] [32] [33], but did not uphold the De Donder 
relation [36], as became clear from discussions with K. 
Sekimoto and D. Lacoste at ESPCI in Paris. It was there 
during a sabbatical leave in 2007 that I first succeeded 
at generalizing BV kinetics to concentrated solutions. I 
spent the next year developing dynamical models of con- 
centrated solutions [35], initially motivated by induced- 
charge electrokinetics [42 . By spring 2009, I synthesized 
the theory in lecture notes [31] and published the BV 
equation for concentrated solutions [35] (Sec. 5.4.2). (See 
also Sekimoto [36].) 

The theory has enabled breakthroughs in modeling 
phase-separating battery electrodes. D. Burch used a 
thermodynamically consistent "Cahn-Hilliard-Reaction" 
model to study intercalation in nanoparticles [13] (with 
a = 0), without connections to the battery voltage. His 
Ph.D. thesis [29 included early simulations of multi- 
particle phase separation and "mosaic instability" [27] 
[28] . Simulations of galvanostatic discharge by P. Bai 
and D. A. Cogswell led to a theory of the suppression 
of phase separation in nanoparticles [15], first reported 
in 2010 [34 , which Cogswell [14 and Stanton [23] ex- 
tended for coherency strain. Meanwhile, T. R. Ferguson 
did the first simulations of phase separation in porous 
electrodes [30] . 

We proceed to summarize the theory and some key pre- 
dictions. First, a thermodynamic framework for reaction 



kinetics is applied to Faradaic reactions, thereby gener- 
alizing the Butler- Volmer and Marcus models to concen- 
trated solutions and solids. The theory is then unified 
with phase-field models by formulating electrode kinet- 
ics in terms of thermodynamic variational principles. Fi- 
nally, the theory is applied to intercalation dynamics in 
LFP nanoparticles. 

II. BUTLER- VOLMER KINETICS IN 
CONCENTRATED SOLUTIONS 

The theory begins with the diffusional chemical poten- 
tial of species z, 

fii = k B T\nCi + /if = k B T]nai (3) 

where q is the concentration, is the absolute chemical 
activity, \l\ x = k B T In ji is the excess chemical potential 
in a concentrated solution, and ji is the activity coeffi- 
cient (di = jiCi). In linear irreversible thermodynamics 
(LIT) [38 , 39 , the flux of species i is proportional to its 
chemical potential gradient, 

F% = -MidVfJn = - A (vci + V^lI) = -D? em Vc, 

(4) 

where Di is the tracer diffusivity, Mi = Di/k B T is the 
mobility (using the Einstein relation), and — V/if* is the 
thermodynamic driving force beyond ideal diffusion (first 
term). By including the mean electrostatic energy (z^ec/) 
below) in the e/ectro-chemical potential, this approach 
has been widely used to describe concentrated solutions 
in electrochemical systems [2 and electrokinetic phenom- 
ena [55] . 

The total flux in [4] has two contributions: (i) free diffu- 
sion driven by concentration gradients (Fick's law), and 
(ii) migration induced by gradients in excess chemical 
potential. In a thermodynamically consistent formula- 
tion of reaction kinetics [31] [36], therefore, the reaction 
complex explores a landscape of excess chemical potential 
\i ex (x) between local minima /if x and ja^ with transitions 
over an activation barrier ji^ x . For the general chemical 
reaction, Si = J2 i SiMi — » ^ SjMj = £2, the activities, 
di = fli a T an d CL2 = Ylj a- 3 , are equal in equilibrium, 
and the forward and backward reactions obey detailed 
balance (R = 0). Out of equilibrium, assuming long- 
lived states with rare transitions (/i| x — \i\ x ^> k B T), the 
net reaction rate is given by 

R = R\^2 — R2-^i 

= V \ Cl e-^t X -^ X )/k B T _ C2e -{n e t x -nT)/kBT 

= V \ e -(v e t x -Vi)/k B T _ e -(^-^2)/k B T^ ^ 

_ v{a\ - a 2 ) 

7t 

which automatically satisfies the De Donder relation [36 , 
Mi — = k B T In (^1-^2/^2-^1 )• Extending Kramer's 
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escape formula [37], the frequency prefactor v depends 
on generalized force constants, ki = d Ji 2 , at the saddle 



dxf 



point and in one minimum (e.g. state 1, with a suitable 
shift of /if*). 

For the general electrode reaction, 



Si = O z 



ne 



R z 



(6) 



we write /ii = \±o + n/i e and \±2 = /i^, where charge 
conservation implies zo — n = zr. We decompose the 
e/ectro-chemical potential into chemical and electrostatic 
contributions, 

Hi = k B T In a,i + Ziefc = (k B T In c* + ziefa) + /if* (7) 

where the term in parentheses describes an ion of charge 
Zie in a mean potential fa in a dilute solution, and 
\i\ x = ksTlnji is the excess chemical potential. The 
electrostatic potential of electrons in the electrode is e , 
while that of the electrolytic solution is <fi. The difference 
is the interfacial voltage, A(f> = <p e — <j). 

Let us consider the broad class of Faradaic reactions 
where the oxidized state is in the electrolyte (e.g .solvated 
Li + ), and the reduced state is either there, or inserted in 
the electrode as a neutral species (e.g. a Li + /e~ polaron). 
The electrochemical potentials of the species in [6] are 



Mo = k B T In ao -\- zoe<p 
H R = k B T\na R + z R e(j) 
H e k B T\na e - efa 



(8) 
(9) 
(10) 



where /i e is the chemical potential of electrons, or Fermi 
level. The electron activity a e = j e c e accounts for the 
concentration of free electrons c e (density of states), as 
well as shifts in the Fermi level via 7 e , due to applied 
stress, phase transformations, donor or acceptor ions, etc. 
as illustrated below for semiconducting electrodes. 

In equilibrium (/ii = ^2), the interfacial voltage, A<p — 
<j) e — </>, satisfies the Nernst equation 



Acj) e 



k B T k a aP e 



(11) 



where n c = n. Out of equilibrium, the current / = neAR 
(where A = electrode area) is related to the local (acti- 
vation, or surface) over-potential, 



r] = Acj) - A(j) e 



(12) 



Specific models of charge transfer correspond to different 
choices of /if*. 

In order to generalize Butler- Volmer kinetics, we ap- 
proximate the mean electrostatic energy of the transition 
state as a linear combination of those of the initial and 
final states: 



state at xr < x$ < xo, then the transfer coefficients, 
a a = 1 — a and a c = a, are related to the symmetry 
factor, 



a 



X% Xr 
Xq ~ Xr 



1 



OLa 



(14) 



which acts as a Br0nsted coefficient [3]. 

Substituting pTp4] into |5j we obtain the BV equation 
□ with 



Io = vneA(coc\ 



,n\l-a a I 
e) C R < \ 



7t 



(15) 



The first factor is the classical formula [TJ [2] , but this is 
only valid for dilute solutions (70 = Jr = 7$ = 7 e = 1)- 
The second factor in braces modifies the exchange current 
for concentrated solutions. In summary, the Faradaic 
current / is related to the interfacial voltage A<p and the 
various chemical activities ao, clr, a e , an d 7j by [T| |TT[ 
[121 and[T5j 



III. MARCUS THEORY FOR 
CONCENTRATED SOLUTIONS 

Just as we have reformulated Butler- Volmer kinetics 
for concentrated solutions, so it is also possible to adapt 
the microscopic theory of charge transfer. The basic idea 
([T]) is that the Faradaic reaction [6] occurs when the excess 
chemical potential of the reduced state, deformed along 
the reaction coordinate by statistical fluctuations, equals 
that of the oxidized state (plus n electrons in the elec- 
trode) at the same point. More precisely, charge trans- 
fer occurs randomly at somewhat lower energies due to 
quantum tunneling [TJ [3] . 

As a first approximation [TJ |3] , we postulate harmonic 
restoring potentials for structural relaxation (e.g. shed- 
ding of the solvation shell from a liquid, or extraction 
from a solid) along the reaction coordinate x from the 
oxidized state at xq to the reduced state at xr\ 



li\ x (x) 

vT(x) 



k B T Indole ) + z oe4> ~ riefa 
+ — (x-x ) 

k B T\n^ R + z R e(j)+ - x R ) 2 (17) 



(16) 



The Nernst equation [TT] follows by equating the to- 
tal chemical potentials at the local minima, fii(xo) = 
I^2{xr) in equilibrium. The free energy barrier is set by 
the intersection of the excess chemical potential curves 
for the initial and final states, /if^ = n\ x (x^) = ^(x^), 
which is an implicit equation for the barrier position, 
x = x$. Using 11 and 12 , we obtain 



ex 

/i t 



k B T\nj t + a a (z e(i) - nefa) + a c z R e</) (13) ner l ~ k B Tln 



cr 
c Q c n e 



If, for consistency, the electric field is constant across the 
reaction coordinate x from xr to xq past the transition 



= T (*t 

= 1^2 {Xr) 

= AG ex 



xof ~ y( x $ 
^T(x ) 



XrY 



(18) 
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FIG. 1: Above: The Faradaic reaction O + ne~ —> R in 
concentrated electrolyte and/or electrode solutions. The 
reactants fluctuate in a landscape of excess chemical 

potential ji ex , shown along a favored reaction 
coordinate x. Charge transfer occurs whenever the 
oxidized state O with n free elections or the reduced 
state R reaches the transition state. Below: Illustration 
of these states for ion reduction by intercalation into a 
solid-solution electrode from a liquid electrolyte. 



where AG ex is the excess free energy change per reaction. 

The overpotential is the total free energy change per 
charge transferred, 



net] 
and 



neAcj) — ksTln 



Q>O0>e 
CiR 



M2O0) - Vi(xr) (19) 



AG = ner\ = AG ex + k B T\n 



cr 
c c n e 



(20) 



is the thermodynamic driving force for the reaction. In 
classical Marcus theory [TJ [4], the overpotential is de- 
fined by ner] = AG ex without the concentration factors 
required by non-equilibrium thermodynamics, which is 
valid for charge-transfer reactions in bulk phases (A~ + 
B — >• A + B~) because the initial and final concentra- 
tions are the same, and thus AG = AG ex = AG (stan- 
dard free energy of reaction). For Faradaic reactions at 
interfaces, however, the concentrations of reactions and 



products are different, and 20 must be used. The missing 
"Nernst concentration term" in |20] has also been noted 
by Kuznetsov and Ulstrup [3] (p. 219). 

Following Marcus, we solve 18 for x$ and relate to 
77. For symmetric relaxation, k^o = kn = we obtain 
the barriers for the cathodic and anodic reactions, 



AG ex = iif 
AG e a x = Lif 



Ll\ X (x ) = 



AG ex 
X 

AG ex 
X 



(21) 
(22) 



respectively, in terms of the reorganization energy, A 



\{x -x R ) 2 . 



by 



Substituting into|5j the current I = neAR is given 



Ioe (ne V f/ik B TX ^ 



-aner} / ksT (l — a)ner)/kB 



(23) 

which has a stronger nonlinear dependence on overpoten- 
tial than in the Butler- Volmer equation. The exchange 
current, 



J = neAe- x l ikBT (c c n e f- a c% 
and symmetry factor, 

1 A fc B T N 



a = 



A 



(24) 



(25) 



are both related to the reorganization energy A. In the 
typical case, A ^> k B T (for thermally stable solvation), 



23 



is 



the symmetric factor is a « ^, and the current 
well approximated by the Butler- Volmer equation up - to 

fairly large overpotentials, |r/| <C \Jl^T m Comparing 
[24] with [l5j we can relate the reorganization energy to 
the activity coefficients defined above, 



^ = k B T\n- ^_ 

4 (707?) 1 " 



IR 



= (l-a)AG eXj raAG ex (26) 
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which is an average barrier for the cathodic and anodic 
reactions, weighted by the symmetry factor. 

In a dilute solution, the reorganization energy Ao 
can be estimated by the classical Marcus approxima- 
tion, Ao = Xi + A , where is the "inner" or short-range 
contribution from structural relaxation (sum over nor- 
mal modes) and A is the "outer", long-range contri- 
bution from the Born energy of solvent dielectric relax- 
ation [US]. For polar solvents at room temperature, the 
large Born energy, A G > 0.5n 2 eV « 20n 2 /c#T, implies 
that single-electron (n — 1), symmetric (a ~ |) charge 
transfer is favored. Many quantum mechanical expres- 
sions for Ao are also available [3]. 

In a concentrated solution, we can use various models 
below to estimate the thermodynamic correction, 7|, due 
to entropic and enthalpic effects in the transition state. 
The exchange current can then be modeled by [T5| with 



the Li + solvation shell) and has reorganization energy Aq, 



7t =7* 



c X /4k B T 



(27) 



where a given by [25] with A = Ao- For example in surface 
adsorption, if the transition state excludes s sites, 7^ = 
(1 — c) _s , then 27 reduces to 



E t /k B T 



(1-c)- 



26 



(28) 



This approach pro- 



where E$ = consistent with 
vides a thermodynamically consistent way to model the 
combined effects of solvation, structural relaxation, en- 
thalpy and entropy on charge transfer kinetics in a con- 
centrated solution, ionic liquid or solid. 



IV. ION INTERCALATION KINETICS 



Let us apply the general theory above to lithium in- 
sertion in a homogeneous solid solution, Li + +e~ — >>Li 



(See 43 below for extensions to heterogeneous solid so- 
lutions with mechanical deformation.) Although lithium 
ions may remain charged, they typically exist as neu- 
tral quasiparticles (Li + M~ polarons [3]) coupled to elec- 
trons localized on solid atoms (M). For example, the stan- 
dard model for LiFeP0 4 [T31IT51 [T8l [41] is a regular so- 
lution [39] [43] for the reduced state R =Li with mole 
fraction c, 



Pl% = k B T\n 



1 



0(1 -2c) 



(29) 



where the first term is the configurational entropy and 
the second the enthalpy of mixing (setting \i = at half 



filling as a reference). The Nernst potential (11) is 



Acf)(c) = Acf) 



MLi(c) , k B T 



\na e (c) (30) 



and, for insertion at the cathode, the open circuit volt- 
age of the battery is Vo(c) = A0(c) — A0 anoc ^ e . If the 
transition state excludes s sites (where s > 1 accounts for 



then using 28 the exchange current 15 takes the form, 



/o(c) = eA(a + a e (c)y-«c«(l-cy-« exp V 

B (31) 

where a + is the Li + activity in the electrolyte and a e (c) 
is the electron activity, which can depend on c, e.g. if Li 
acts as a donor impurity in a semiconducting insertion 
electrode. 

This example illustrates how the theory unifies charge- 
transfer kinetics with battery thermodynamics. In 
porous electrode models for Li- ion batteries [2] [5] [7] , the 
rate expression [2] is assumed, which corresponds to [3l] 
with a + = a e = 1, O = 0, and s = 1. The dilute 
electrolyte activity, a + = c + has also been included [7], 
but these models are inconsistent with the use of con- 
centrated solution models (a + 7^ c + ) to fit experimental 
data. A more serious inconsistency is the fitting of the 
open circuit voltage, Vb(c), which should be related to 
the local ionic activities [29] in both the Nernst equation, 
[30] and the exchange current density, [3l] It is also in- 
consistent to model solid diffusion by Fick's law with a 
fitted nonlinear diffusivity, rather than using ionic activ- 
ity gradients in [4] 



V. VARIATIONAL FORMULATION OF 
ELECTRODE KINETICS 

In order to describe phase separation, we adopt the 
Cahn-Hilliard (CH) model, based on the Helmholtz free 
energy functional, 



F[c] 



Psf(c) 



1 



Vc • kVc 



dV 



(32) 



where n is the "gradient penalty tensor" [39, 43]. (We 
define / and c per crystal site, and introduce the site 
density p s .) The bulk chemical potential is defined vari- 
ationally, 

1 5F V-^Vc 

V = —jt = f ( c ) ( 33 ) 

Ps 0C p s 

as the free-energy change for adding a bulk "particle" 
(delta function 5c). If f(c) has multiple minima rep- 
resenting distinct stable phases, then the gradient term 
counteracts the tendency to phase separate and leads to 
a diffuse interface, whose equilibrium profile satisfies \i = 
constant. 

There were several breakthroughs in 2004 that set the 
stage for our work. Han et al. [41] first applied the CH 
equation to model diffusion in LFP using the regular solu- 
tion model (29). Garcia et al. [44 formulated variational 
principles for electromagnet ically active systems, which 
unify the CH equation with Maxwell's equations. Guyer 
et al. [45 first represented a diffuse electrode/electrolyte 
interface with a continuous phase field. They described 
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the kinetics of electrodeposition [46] by linear Allen-Cahn 
kinetics [39], / ex A/i, but did not make connections with 
charge-transfer theory. 

Our theory followed from these developments and the 
observation of phase boundaries at LFP surfaces [40j [47] . 
The exponential (Tafel) dependence of the current on the 
overpotential, defined in terms of A/i in [33] was first re- 
ported in 2007 jTSJ H2J ESi , but with a spurious pre-factor, 
noted by Burch [13 , 29]. The complete theory, connecting 
the CH model to the BV equation and the battery volt- 
age, appeared in 2009 [31] and was later applied to LFP 
nanoparticles [14} [15]. For lithium intercalation, the re- 
duced state has zero charge (zb,= 0), so the original CH 
functional 



Lf intercalation 



32 



with /(c) given by [29] was used to construct 
the electrochemical potential [33] in the Butler- Volmer and 
Nernst equations. 

For the general Faradaic reaction, 

O + ne~ = s i°T + ne ~ Yl s i R 7 = R ' ( 34 ) 



we define all the chemical activities and diffu- 
sional electrochemical potentials pk variationally from 
the Gibbs free energy G: 



Po = k B T\na + zoe^o ^ 

i 

Pr = k B T\na R + z R e<j> R = ^ 



/c#Tlna e 



j 

1 SG 

Ps $Ce 



Si SG 
Ps SCi 

Sj SG 

Ps SCj 



(35) 
(36) 
(37) 



The Nernst voltage A(f) eq is defined variationally by set- 
ting jjL R = jio + np e . According to [TT] [12] and [20] the 
activation overpotential 



Empty 
Channels 



Intercalated 
/Channels 




FIG. 2: Insertion kinetics in reaction-limited LFP 
nanoparticles. (a) The low-energy phase boundary 
between FeP04 and LiFeP04 propagates as an 
"intercalation wave" [TT] El [HI US [33 across the {010} 
facet, filling the crystal layer by layer, as the reaction 
rate (43) is enhanced by concentration gradients and 



elastic stress. [ Bai et al. [15 .] (b) Collision of 
intercalation waves nucleated by lithium wetting of the 
side facets prior to complete filling. [D. A. Cogswell, 
unpublished] 



V 



PR-(po+riPe) P2-P1 



ne 



ne 



(38) 



is the change in total chemical potential of the Faradaic 
reaction per charge transferred (n c = n). In this way, the 
diffusional electrochemical potentials pi = s f c enter the 
Butler- Volmer ([!]) and Marcus (23) equations, as well as 



the exchange current 15 



The following free energy captures the key physics of 
ionic materials [IMS] M M S3 HE] : 



G = I [Ps{m+Pe(i)) 
1 



(Vc • ftVc - Vcj) ■ s p V(j) + a : e) 
7s(n,c, 4>)dS 



dV 



(39) 



where c is the set of concentrations, e p the permittiv- 
ity tensor, k the gradient penalty tensor, a the stress 
tensor, e the strain tensor, 7 S the surface tension, and 
p e = ZiCCi the bulk charge density. The electrostatic 



potential acts as a Lagrange multiplier controlling the to- 
tal ion densities [14] |44] while enforcing Maxwell's equa- 
tions [48 . Assuming linear elasticity, the stress is given 
by Hooke's law, = Cijki^kU where C is the elastic con- 
stant tensor. The lattice-preserving "coherency strain" , 



1 / dui_ 

2 1 dx j 



3uj_ 
dxi 



ijm c m 



(40) 



is the total strain due to compositional inhomogeneity 
(first term) relative to the stress-free inelastic strain (sec- 
ond term). In a mean- field approximation (Vegard's law), 
each molecule of species m exerts an independent strain 
(lattice misfit between c m = 0, 1). 
Through the electro-chemical potentials, 



Pt 



dci 



k B T\nai 



(41) 

the Nernst voltage [TT] and exchange current [15] have 
novel dependencies on concentration gradients and elas- 
tic stress, which lead new modes of ion intercalation (J2|. 
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For ion insertion 0, we can model the transition state 

1 / A r , (7 : ej 



as 



7t = (1 - °i) * ex P 



Ao 
4 



(42) 



where is the transition-state strain, s > 1 the number 
of excluded sites, and Ao the reorganization energy. For 
insertion in a regular solid solution, the exchange current 
is 

io = /o(c)exp^ (43 ) 

where 7p(c ) is the limit of zero stress and uniform con- 
and = — ae°, is the "activation 



centration 31 
strain" [49]. 

From the LIT flux [4] conservation of mass leads to the 
CH equation, 



Ri 



where F 7 - 



■^V^. (44) 
k B Tp s Sci 



The homogeneous reaction rate Ri, describing bulk 
charge transfer reactions (e.g. water self-ionization) , 
is often neglected, but could be described by our the- 
ory. The current density is also expressed variationally, 
J = J7 ZieFi. The CH equation is fourth-order and re- 
quires two boundary conditions at each interface. One 
relates the normal flux to the reaction rate, h • = Ri, 
where Ri = for an inactive species and Ri = for 
an active species. The other is the "variational boundary 
condition" , 



[n • nVci) = 



dci 



(45) 



(where [X] denotes the jump in X at an interface), which 
ensures continuity of the chemical potential [13] and con- 
trols surface wetting and nucleation [T5] . 



VI. ION INTERCALATION WITH PHASE 
SEPARATION 

Size-dependent solubility. It is well known that 
the spinodal gap of unstable compositions in a solid so- 
lution shrinks with decreasing particle size [50], as low 
energy long- wavelength fluctuations are prohibited. The 
miscibility gap also shrinks, since the relative cost of 
a phase boundary increases with the surface to volume 
ratio [5j. Elastic coherency strain enhances both ef- 
fects by increasing the energy penalty for phase separa- 
tion nana ed. 

The theory predicts novel effects of kinetics on solu- 
bility, due to ion exchange with the electrolyte reser- 
voir p~3j[14]. At constant voltage, the spinodal and misci- 
bility gaps shrink with increasing exchange current [13], 
since fluctuations with fast reactions promote unfavor- 
able concentration gradients. Experimental solubility 



data [17 can be interpreted using a ID Cahn-Hilliard- 
reaction (CHR) model Q~3j[29], or a more accurate, 2D 
reaction-limited model with coherency strain [T4] . 

Intercalation waves. In anisotropic crystals, ion 
intercalation dynamics can be limited by diffusion or re- 
actions, independently in each direction [18 . In LFP 
nanoparticles, the concentration tends to remain uniform 
in the [010] direction due to the fast diffusion [19] [52] 
and elastically unfavorable phase separation [14] in short 
one-dimensional channels, unblocked by Fe anti-site de- 
fects [20J. After [010] depth- averaging [18 , the full CHR 
model reduces to a new, reaction-limited model on the 
active {010} crystal facet, 



(46) 



(where Jo = Io/A s is the current per surface site), which 
is an electrochemical generalization of the Allen-Cahn 
equation [39 (EAC) with with a nonlinear dependence 
on the variational overpotential, 



rj 



ep s oc 



(47) 



where f(rf) = I/Io can come from the BV ([I]) or Mar- 
cus ( [23] ) equations. For lithium insertion in a regular 
solution [14][T5], the exchange current Iq(c, V • ftVc, a) is 



given by 43 



The theory has led to a paradigm shift in understand- 
ing insertion dynamics. In contrast to prior battery mod- 
els, which assume a uniform reaction rate over the par- 
ticle surface, limited by radial solid diffusion [2] [5] or 
a shrinking-core phase boundary [6j [7], the EAC limit 
of the CHR model predicts reaction-limited intercala- 
tion waves (or "domino cascades" [53 J, sweeping across 
the active facet, filling the crystal layer by layer ([2| 
[HJ[T5j[l8> 33, 54 . Intercalation waves result from nu- 
cleation at surfaces or defects, or by spinodal decom- 
position [15], and trace out the voltage plateau at low 
current, as shown in [3] Experimental evidence of one- 
dimensional nucleation and growth, consistent with the 
EAC theory [15] [18] , has recently been reported [26] . 

Driven phase separation. The theory has led 
to some surprising conclusions about electrochemically 
driven phase transformations in nanoparticles. Our first 
paper included a proof that traveling-wave solutions of 
the EAC equation only exist over a finite range of thermo- 
dynamic driven force p~8j [32] , but the implications were 
unclear without the connection to the battery voltage. 
Once the generalization of B V kinetics was achieved [31] , 
simulations and analysis revealed how phase separation 
is strongly influenced by the applied current [14] [15] . As 
suggested by Malik et al. [16 based on bulk free energy 
calculations, phase separation is suppressed during bat- 
tery operation, but the mechanism is related to activation 
overpotential at the surface. 

As explained by Bai et al. [15], phase separation is sup- 
pressed as the applied current approaches the exchange 
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FIG. 3: Galvanostatic discharge into a single LFP 
cathode nanoparticle without coherency strain. Top: 
Linear stability (spinodal) diagram for the homogeneous 
state versus current / (scaled to the exchange current at 
X = 0.5) and mean Li filling X. (a) Battery voltage 
versus state of charge (X) at different currents /. (b) 
Evolution of the concentration profile at low current 
I = 0.01 with insertion waves nucleated at the wetted 
[001] side facets, (c) Quasi-solid solution behavior close 
to the critical current / = 0.5, showing homogeneous 
filling ahead of the waves. [ Bai et al. [T5] ] 



current. A linear stability analysis of the EAC model 
for BV kinetics with a = \ shows that concentration 
fluctuations of wave number k grow at a rate [15] 



s(k- c, I) = -\l J + ^r(fj f + Kk 2 ) + / 



Jo 



1 



43 



(48) 

where 77(c) and Jo(c) are the activation overpotential and 
exchange current density, respectively, in a homogeneous 
state of concentration c sustaining a current /. The range 
of linear instability, s(c, /) > 0, predicted by 48 and 
with elastic strain is shown in[3ja). 

The theory predicts a critical current, of order the ex- 
change current, above which phase separation by spin- 
odal decomposition is impossible. Below this current, 
the homogeneous state is unstable over a range of con- 
centrations (smaller than the zero-current spinodal gap), 
but for large currents, the time spent in this region is too 
small for complete phase separation. Instead, the parti- 
cle passes through a transient quasi-solid solution state, 
where its voltage and concentration profile resemble those 
of a homogeneous solid solution. 

The suppression of phase separation predicted by the 
theory is very general. When nucleation is possible (see 
below), a similar current dependence is observed. As 
shown in [3j the bulk crystal fills ahead of the surface- 
nucleated intercalation waves in the quasi-solid solution 
regime, while homogeneous filling occurs above the crit- 
ical current. With elastic coherency strain included, the 
suppression of phase separation is analogous, but even 
stronger [14]. 

In the case of LFP, the critical current is an order of 
magnitude smaller than the exchange current [14]. The 
surprising implication is that phase separation does not 
occur in situ during normal battery operation, which 
helps to explain the high-rate capability and extended 
lifetime of LFP nanoparticles, due to larger active area 
and smaller elastic stresses. On the other hand, phase 
separation occurs at very low currents and can be ob- 
served in particles trapped at intermediate compositions 
for ex situ measurements Q . 

Phase separation patterns. For quantitative in- 
terpretation of experiments [14], it is essential to ac- 
count for the elastic energy. Even a simple elastic model 
can dramatically improve the inference of lithium dif- 
fusivity from intermitted titration with phase separa- 
tion [56 . Crystal anisotropy generally leads to striped 
patterns in equilibrium [2TH23] , The stripe spacing fol- 
lows from the balance between elastic energy (favoring 
short wavelengths at a stress-free boundary) and interfa- 
cial energy, favoring long wavelengths to minimize inter- 
facial area [13] . In layered crystals, the lattice misfit can 
cause tension and compression in different directions, and 
Stanton and Bazant [23] predicted that this could cause 
phase boundaries to tilt with respect to the crystal axes. 
For example, lithium insertion in LFP causes contraction 
of the iron phosphate planes in the [001] direction and ex- 
pansion in the [100] and [010] directions [40]. Cogswell 



FIG. 4: Simulations of reaction-limited phase separation of a 500nm single crystal of Lio.sFePC^ into Li-rich (black) 
and Li-poor phases (white) in an electrolyte bath at zero current and zero pressure, consistent with ex situ 
experiments [40j [471 ES] • ( a ) Coherent phase separation, where lithium insertion causes contraction along [001] (c) 
axis and expansion along the [100] (a) axis, leading to tilted interfaces aligned with {101} planes, (b) Loss of [001] 
coherency (e.g. due to microcracks [40]) causes the phase boundaries to rotate to align with {100} planes. 

[ Cogswell and Bazant [14].] 



and Bazant [14] explained the experimentally observed 
phase morphologies [47] [55] and used the theory to in- 
fer the gradient penalty n from the stripe spacing. As 
shown in Fig. [4] the model enables realistic simulations 
of reaction-limited phase separation. 

Surface wetting and nucleation. The theory also 
predicts how surface thermodynamics affects reaction ki- 
netics. Using the boundary condition [45] for the EAC 
equation [46] the nucleation of intercalation waves in LFP 
nanoparticles due to wetting of the {100} facet can be 
simulated with ^p) or without (J^) coherency strain [15] . 
The nucleation barrier (critical voltage past the equilib- 
rium plateau to observe a sudden current) can be pre- 
dicted as a function of particle size, in excellent agree- 
ment with experiments (D. A. Cogswell, unpublished). 
The effect of surface wetting on intercalation kinetics may 
be responsible for the surprising performance boost from 
phosphate-glass thin films on active nanoparticles of LFP 
[10] and Li^CoOs [57]. 

Mosaic Instability. The complex dynamics of ion 
intercalation in nanoparticles also affects the macrosopic 
response of a composite electrode. Using the CHR 
model [13] for a collection of particles in a reservoir of 
constant chemical potential, Burch [29] discovered what 
we called the mosaic instability, whereby particles filling 
uniformly at constant total current suddenly fill one-by- 
one. The mechanism is simple: As the particles enter 
the miscibility gap, the first to phase separate will sud- 



denly lower its chemical potential and sustain a much 
larger current until it becomes full, at which point the 
current is shared by the others until the process repeats. 
Interestingly, in the CHR model, the smaller particles 
transform first, because diffusion brings the core concen- 
tration to the spinodal more quickly. In contrast, for 
reaction-limited particles, the larger ones are more likely 
to nucleate first. 

Meanwhile, Dreyer et al. [27 described the mosaic 
instability and demonstrated it experimentally in two 
clever ways. They connected it to the zero-current volt- 
age gap between charge/discharge cycles in LFP batter- 
ies, as well as the dynamics of an array of balloons con- 
nected to a single pressure line. Their theory [27] [28] 
is similar to the "pseudo capacitor limit" of our models, 
where each particle maintains uniform (but time depen- 
dent) concentration, without internal phase separation. 

Porous electrodes. The theory has recently been ex- 
tended to porous electrodes undergoing phase transfor- 
mations [30 using classical porous electrode models to 
describe the electrolyte [2]. Diffusion in the electrolyte 
provides a spatiotemporal bias for the mosaic instabil- 
ity, so that narrow fronts of stochastic particle phase 
transformations propagate away from the separator (dis- 
charge) or current collector (charge) at low current 
The width of the reaction front separating low-density 
and high-density phases increases with current until the 
entire electrode becomes active, which is a macroscopic 
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(a) (b) slow C/50 charge c s {x,y,t) 




Filling Fraction x/i_ x x/i_ x 

FIG. 5: Finite- volume simulations of macroscopic phase separation in a porous LFP cathode [30]. (a) Voltage versus 
state of charge for slow charge/discharge and fast discharge, (b)-(d) Profiles of mean Li concentration in the solid at 
the red points in (a), showing narrow reaction fronts with mosaic instabilities, associated with voltage fluctuations, 
(e) Electrolyte Li concentration for fast discharge (d), showing diffusion limitation. [Ferguson and Bazant [30] ] 



analog of the suppression of phase separation in nanopar- 
ticles. 

The theory can predict experimental data for phase- 
separating porous electrodes for the first time without 
artificially fitting voltage plateaus or enforcing phase 
boundaries. Using the EAC regular-solution model for 
LFP [14] [15] in a porous electrode [30], the model can 
predict tiny voltage oscillations (< /c#T/e =25mV) and 
the limiting "voltage gap" at very low currents observed 
in experiments [27], which give way to a smoothly de- 
caying voltage at high current [10] Nucleation is 
also important at the porous electrode scale [26] and can 
be naturally included in theory. Using a two-component 
regular solution model for graphite intercalation [31] , the 
theory is also able to fit experimental data for diffusion- 
limited intercalation of long graphite anodes [58] with 
multiple reaction fronts of different colors corresponding 
to different stages of lithiation. 

VII. CONCLUSION 

This Account describes the development of a theory of 
electrochemical kinetics in condensed matter. The theory 
has shed light on intercalation dynamics in advanced Li- 



ion batteries with phase-separating nanoparticles. Other 
possible applications include oxygen electrocatalysis in 
fuel cells and electrochemical stability of ionic liquids. 
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